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We study the order-from-disorder transition and reentrant magnetism in La2-iSr^Cui_ z Zn z 04 
within the framework of a long-wavelength nonlinear sigma model that properly incorporates the 
Dzyaloshinskii-Moriya and XY anisotropies. Doping with nonmagnetic impurities, such as Zn, is 
considered according to classical percolation theory, whereas the effect of Sr, which introduces charge 
carriers into the Cu02 planes, is described as a dipolar frustration of the antiferromagnetic order. 
We calculate several magnetic, thermodynamic, and spectral properties of the system, such as the 
antiferromagnetic order parameter, <7o, the Neel temperature, Tjv, the spin-stiffness, p 3 , and the 
anisotropy gaps, Adm and Axy, as well as their evolution with both Zn and Sr doping. We explain 
the nonmonotonic and reentrant behavior experimentally observed for Tm(x,z) by Hiicker et al. in 
Phys. Rev. B 59, R725 (1999), as resulting from the reduction, due to the nonmagnetic impurities, 
of the dipolar frustration induced by the charge carriers (order-from-disorder). Furthermore, we 
find a similar nonmonotonic and reentrant behavior for all the other observables studied. Most 
remarkably, our results show that while for x « 2% and z = the Dzyaloshinskii-Moriya gap 
Adm = 0, for z — 15% it is approximately Adm ~ 7.5 cm -1 . The later is larger than the lowest 
low-frequency cutoff for Raman spectroscopy (~ 5 cm -1 ), and could thus be observed in one-magnon 
Raman scattering. 

PACS numbers: 



I. INTRODUCTION 

High temperature superconductivity is obtained by in- 
troducing charge carriers (holes or electrons) into Cu- 
based Mott-Hubbard insulators. At half filling such sys- 
tems exhibit long range antiferromagnetic (AF) order in 
the ground state, which is however rapidly suppressed 
by the introduction of charge carriers. Consider, for ex- 
ample, the case of La2Cu04, the simplest of the par- 
ent compounds. The replacement of La 3+ by Sr 2+ ions 
in La2_ £C Sr :c Cu04, through which holes are doped into 
the CuC-2 planes, causes the destruction of the canted 
Neel order already at x 0.02ji showing that the doped 
holes strongly frustrate the underlying AF order within 
the CuC-2 planes. Conversely, doping of iso-valent non- 
magnetic impurities into La2Cu04 is known to have less 
dramatic effects. For example, by replacing Cu 2+ for 
Zn 2+ in La2Cui_ z Zn z 04fSi the AF order is suppressed at 
much higher impurity concentration and the monotonic 
and rather smooth decrease of the Neel temperature Tjv 
with z can be described, at least in the limit of low dilu- 
tion, within classical percolation theory. 

The combined effect of adding holes and nonmagnetic 
impurities into L^-ajSrzCui-zZn.zC^ was studied exper- 
imentally by Hiicker et alfi It has been found that, even 
though each kind of impurity independently suppresses 
the AF order, an enhancement of the antiferromagnetism 
can occur when these impurities are combined. For ex- 
ample, in the presence of Zn, z > 0, the AF order survives 
for x > 0.02.8 More interestingly, it was also found that 
for x — 0.017, the Neel temperature exhibits a nonmono- 
tonic behavior as a function of Zn doping. In fact, TV is 



first enhanced from 125 K at z — to 144 K at z = 0.05, 
and only then starts to be suppressed by the dilution (see 
Fig. 4 of Ref.0). 

This remarkable reentrant effect was immediately ad- 
dressed theoretically by Korenblit et al.fi who first sug- 
gested that Zn reduces the frustration induced by Sr, 
through the dilution of ferromagnetic bonds. These au- 
thors considered an isotropic Heisenberg spin system and 
neglected the Dzyaloshinskii-Moriya (DM) and pseudo- 
dipolar (XY) anisotropies, including only the much 
smaller interplanar superexchange J±. More recently, 
however, the same authors suggested that J± has little 
or nothing to do with the suppression of the AF order 
in La2-xSr x Cui_ z Zn z 04 ( J± was shown to be almost x 
and z independent), which should be, instead, solely de- 
termined by intraplanar correlations. 5 This is in agree- 
ment with the recent findings of Juricic et al.fi where it 
has been shown that the robustness of the canted Neel 
state is determined by the DM gap. The collinear long 
range AF order is destroyed when the renormalized DM 
gap vanishes, at x w Adm/ J ~ 0.02, and for higher 
doping the magnetism becomes incommensuratei&i The 
DM and XY anisotropies have also been shown to be 
behind the unusual magnetic susceptibility response in 
La2- a; Sr a ;Cu04, s - 9 for a rather wide range in doping and 
temperature. In addition they are also responsible for the 
appearance of a field-induced mode in the one-magnon 
Raman spectrum of L^-gSrsCuCU i 1 ^ 11 Thus it is clear 
that any realistic description of the reentrant magnetism 
in La2- a; Sr a ;Cui_ z Zn z 04 must properly take into account 
such anisotropies. 

In this paper we revisit the problem of the reduction 
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of frustration in La2- 3; Sr 3 ;Cui_ 2 Zn z 04, first proposed 
by Korenblit et al.fi within the framework of a long- 
wavelength nonlinear sigma model (NLSM) that properly 
includes DM and XY anisotropies, dilution and frustra- 
tion. We show that indeed dilution weakens the frus- 
tration by reducing the dipolar-magnon coupling con- 
stant. The result is a nonmonotonic behavior not only 
for Tjv, but also for other observables like: order param- 
eter, weak-ferromagnetic moment, and anisotropy gaps. 
Moreover, since the destruction of the AF order is di- 
rectly related to the vanishing of the DM gapj& we find 
that the reentrant antiferromagnetism is accompanied by 
a reentrant behavior for the DM gap itself. Finally, we 
show that the NLSM formulation, together with classi- 
cal percolation theory, can also describe the experiments 
in the highly Zn-diluted limit, provided the appropriate 
bond percolation factor is used. 

The paper is organized as follows. In section[H]we dis- 
cuss the different impurities to be considered, we review 
the derivation of the anisotropic NLSM for La2CuC>4, in- 
troduce dilution by Zn according to the classical percola- 
tion theory, and include the dipolar frustration within the 
framework of the Shraiman and Siggia modeled In sec- 
tion lllll we present our general model and compute all the 
renormalizations of the physical quantities. In section lTVl 
we discuss the effect of the reduction of frustration in dif- 
ferent physical quantities and propose experiments that 
could lead to the observation of these effects. This section 
is written in such a way that the reader not interested 
in the theoretical details can skip sections ITT1 and ITTTI In 
section[^we present our conclusions. In AnnendixlAlwe 
review the continuum limit for the site and bond perco- 
lation factors, while in Appendix[B]we provide a detailed 
derivation of the NLSM with dilution. 



II. IMPURITIES IN HOST La 2 Cu0 4 

La2CuC>4 is a layered antiferromagnet with a rather 
large XY anisotropy which below 530 K is in the low- 
temperature orthorhombic (LTO) phase. Although the 
orthorhombicity resulting from the staggered tilting of 
the CuC>6 octahedra is small, it is responsible for the pres- 
ence of an antisymmetric exchange of the DM type once 
the spin-orbit coupling is considered. Together, these two 
anisotropies, DM and XY, gap the transverse spin-wave 
excitations along the a and c orthorhombic axis respec- 
tively, and lead to a canted antiferromagnetic ordering 
along the b orthorhombic direction, see Fig.^Jj). 

Theoretical studies of La2CuC>4 concentrate on the Cu 
lattice, because all the other atoms (La 3+ , 2 ~) are in 
a closed shell configuration. In the crystal Cu 2+ has a 
missing electron in the d x 2_ y 2 orbital, which carries a 
spin 1/2 and orders antiferromagnetically below Tn=325 
K. There are two kinds of impurities which, at some 
value of doping, destroy the antiferromagnetic order in 
the layered La2CuC>4 compound. Replacing La 3+ with 
Sr 2+ leads to a rapid suppression of the antiferromag- 
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FIG. 1: Left: Tetragonal (xyz) and orthorhombic (abc) co- 
ordinate systems. Right: In the LTO phase the tilting axis 
of the CuOe octahedra, represented by the vector D + (in 
green), determines the orthorhombic a axis. The DM and 
XY anisotropies then determine that the staggered order pa- 
rameter n (in red) is oriented along the orthorhombic b di- 
rection while the weak-ferromagnetic moment L (in blue) is 
perpendicular to the Cu02 plane. 



netic order in the plane. This happens because the holes 
donated to the planes form Zhang-Rice singlets with the 
local moments of Cu 2+ ji£ which act as effective mobile 
holes in the spin lattice. Doping the crystal with Zn, 
i.e. substituting Cu 2+ with iso-valent Zn 2+ , introduces 
vacancies at the Cu positions. In fact, Cu 2+ has an elec- 
tronic configuration [Ar]3d 9 , while Zn 2+ is [Ar]3d 10 . This 
means that there is no magnetic moment at the Zn po- 
sition and therefore the effect of Zn-doping is simply to 
dilute the spins within the plane. 

The effect of the competition between Sr and Zn im- 
purities in La2- a: Sr ;f Cui_ z Zn z 04 was experimentally in- 
vestigated in Ref. y and addressed theoretically in Refs. 
H and [3 . Both theoretical works have, as starting 
point, an isotropic quantum Heisenberg Hamiltonian on 
a square lattice. However, recent experimentalise and 
theoreticalSiiiiiSiiLiS studies have emphasized the very 
important role of DM and XY anisotropies and thus we 
shall now describe the effect of dilution and frustration 
in an effective model including anisotropies. 



A. Undoped La2Cu04 

At long wavelengths, the isotropic Heisenberg Hamil- 
tonian is well described, in the paramagnetic phase, by 
a 0(3) NLSMii^ In this paper, however, we use a gen- 
eralized NLSM, derived by Chovan and Papanicolaou in 
Ref. and Silva Neto et al. in Ref. H which includes 
the DM and XY anisotropies. The action of the model 
reads 

Stot = J dr[LH + Ldm + Lxy + Lwz], 



<ho> 

iS Sflj ■ (ttj x dottj). 

je2Dlattice 
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In Eq. JQ, Lh, Ldm, Lxy and Lwz are, respectively, 
the Heisenberg, DM, XY, and Wess-Zumino terms, ft 
is a unit vector along the spin direction, S = SSI, and 
S = 1/2. The super-exchange parameter is denoted by 
J. The DM vector and the anisotropy matrix I\j are 
defined on the Cu-Cu bondsi22i2I In the (x, y) coordinate 
system, see Fig. Q] T> x = (0,d,0), T) y = (d,0,0), f x = 
diag(Ti + T 2 , T± - T 2 , T 3 ), and T y = diag(T 1 - T 2 , T x + 
Ta, r 3 ). Here, the index x means along the horizontal Cu- 
O-Cu bonds, y means along the vertical Cu-O-Cu bonds, 
d ~ 10- 2 J, ri )2|3 ~ 10- 4 J, and Ti > T 3 . 

Our next step is to separate the fast 1 and slow n vary- 
ing spin components 

Oi = (-1)^ VI - K) 2 + ah 

~ (-ly^ + aU-^^a 2 !!^. (2) 

In Eq. (|2l o is the lattice constant and i is an index on 
the 2D lattice, i = {p,q); (—1)* should be understood as 
(— l) p+q , rii as n P;g . After substituting Eq. Q into Eq. 
@J we find 

Stot = Jd 2 r{^ [(Vn) 2 + 81 2 ] + • (n x 1) 

+ ^! (ri -r 3 )n 2 -^l.(nxri)}, (3) 
where d + = (D^ + B y )/2. 

B. Effect of non-magnetic impurities - Zn doping 

Recently, an effective field theory for the Heisenberg 
antiferromagnet with non-magnetic impurities was de- 



where we have defined Z — i4T 2 /P 00 , g is the usual NLSM 
coupling consta nt, and c is the spin-wave velocity. Here 
we used that y / 2gcS 2 /Ja 2 d + = 2\f2 Sdel = m a e^ and 
(4gcS 2 /a 2 )(T 1 - T 3 ) = 32JS' 2 (ri - T 3 ) = m 2 , where 

is the unit vector along the a orthorhombic direction. 
Therefore, the last two terms in Eq. (@J correspond, re- 
spectively, to the DM and XY anisotropy gaps, show- 
ing that the spin ordering has an easy-axis along the or- 
thorhombic b direction, see Fig. 

From the above action it is clear that when only 
Zn impurities are doped into La2Cui_ z Zn z 04 the two 
anisotropy gaps renormalize according to 

M QjC = VZ m atC , (5) 



rived by Chen and Castro Netoi 2 ^ While in Ref. ^ the 
authors only considered the isotropic case, here we gen- 
eralize the results of Ref. 112 to include the DM and XY 
anisotropies. 



We introduce a function pi, which is 1 on a Cu site 
and zero on a Zn site, with the property p 2 — Pi. We 
then assume this function to be smooth and expand pj 
in the neighborhood of its nearest neighbor site. In a 
continuum limit, pt and Kij = piPj are replaced, respec- 
tively, by p(r) and K(r), which are the so called site and 
bond percolation factors. As it is shown in the Appendix 
El for the case of a homogeneous distribution of static 
impurities, they can be replaced by their average values, 
Poo(z) = 1 — z and K{z) = 1 — 3z. In the following 
we omit the z-dependence in Poo and K for the sake of 
shortening the notation. 



In the presence of dilution, the action (QJ will be mod- 
ified as follows: the terms L H ,L DM , and L X y, which 
contain J2<ij> /(^i, wm ^ e multiplied by K. The 
Wess-Zumino term, on the other hand, contains J7f and 
therefore it is proportional to p 3 (r) = p(r). Thus, within 
the homogeneous distribution assumption, it will be mul- 
tiplied by Poo. In this way, K and Poo appear simply as 
prefactors of the relevant integrals^ A detailed deriva- 
tion of the NLSM in the presence of nonmagnetic impu- 
rities is given in Appendix [bJ 



The next step is to integrate out 1 in Eq. <PJ to obtain 
the action in the presence of dilution (see Appendix [Bj 



(4) 

I 

and thus decrease rather smoothly with dilution and van- 
ish at the percolation threshold. 

In what follows we shall switch freely between Adm 
and M a , when refering to the DM or in-plane gap, and 
also between Axy and M c , when refering to the XY or 
out-of-plane gap, withouth any loss of generality. 



C. Derivation of the Neel temperature 

Let us now derive an expression for the Neel temper- 
ature in terms of the anisotropy gaps and spin-stiffness, 
renormalized by dilution. Starting from the action in 
Eq. iQJ, we split the n-field into its longitudional ctq 



I 

Szn = 2gc( K /Poo) J dT J d 2 r{(d T n) 2 + Z[c 2 (Vn) 2 + m 2 n 2 + m 2 c n 2 c ]} , 
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and transverse ni components, n = (n a ,(ro,n c ). Here 
= (n a ,n c ), and ao = const, is the order parameter. 
The action then reads 



(6) 



with 

A d -- 



(P oc /K)u; 2 n + Kc 2 k 2 +dmg[Kmi,Kmi], (7) 



where uJ n — 2nn/(3 are the Matsubara frequencies, 1 2 
is a 2D identity matrix, and (3 = l/fc^T is the inverse 
temperature. Defining Tr as ^ w J" ePk/(27r) 2 and 
introducing the fixed length constraint into the action 
through a Lagrange multiplier, A, the partition function 
can be written as 

Z = J Dn S(n 2 - P^) exp{-S* d [n]} 

= / Da Dn ± D\cxp ( - J—Tr[n ± A d n± 
J I- 25c 

= J Da DXe X p{-S ef f[\,ao]}, (8) 
where the effective action reads 

S ef f[X,a ] = ±-Tr[iX^-iXP 00 }+^ ]T Trln[^ Q+l A]. 

a—a.c 

(9) 

In Eqs. ((HI and J2j we used that the average length of 
spin per lattice site is Poo in the diluted case. 22,23 

The thermodynamic properties of the diluted system 
can be determined, at the mean field level, by solving the 
saddle point equations of the effective action From 
5S e ff[X,aa)/SX — we find 



E Tr^T 1 " 

a— a,c 



iA 



(10) 



The only condition for the trace to converge is to choose 
iX = m,Q, a real positive number. As usual, we interpret 
it as the inverse correlation length, uiq = l/£. 

At the Neel critical point the order parameter ao van- 
ishes, the correlation length £ diverges, and the equation 
for the critical temperature is 



gcK 



EE 



a—a.c uj n 



d 2 k 



1 



(2tt) 2 W 2 +Zc 2( k 2 +m 2/ c 2)' 



After performing the Matsubara summation and the in- 
tegration over momenta, the last expression simplifies to 



= 0, (11) 



where Pn = and we defined the renormalized 

spin stiffness for the diluted system 



Ps = c 7T 



1 



1 



2g g c KP a 



(12) 



with g c defined through the ultra-violet momentum cut- 
off A = 47r/<7 c , as usual. 

T M (0,Z) dnK) 
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FIG. 2: Neel temperature (in K) as a function of Zn con- 
centration for x = up to the highly diluted regime. The 
experimental data for either powder and crystal samples was 
taken from Refs. Q and0 

The bond percolation factor K{z) — 1—32, proposed in 
Ref. I22I (see Appendix El describes well the experimen- 
tal data for Tjq(z) only in the low doping rangei 3 ^ 2 ^ For 
highly diluted samples, however, single crystals were not 
available before the work reported in Ref. 0- The exper- 
iments were usually performed with single crystal sam- 
ples below the doping threshold of z ~ 20%, while above 
25% powder samples were used. After Ref.|2J, which used 
both single crystals and powder samples, our conclusion 
is that for the heavily doped samples the bond perco- 
lation factor should be modified according to Refs. I24I 
andl25l which suggest that K(z) = 1 — nz + nz 2 /2. In 
this Watson-Leath (WL) prescription, Tn is zero at the 
percolation threshold, z p . Theoretical, experimental, and 
numerical studies of highly diluted La2Cui_ z [Mn,Zn] z 04 
are presented in Ref. 0- There the authors describe an 
experiment where the critical point is found at z = 42%. 
As it is observed in Fig. using the bond percolation 
factor as in the WL prescription, we obtain the correct 
description in the heavily doped regime. 

For the plot in Fig. [2 we used J = 100 meV, which is 
much smaller than the real value of s» 135 meV, but is 
the value that gives T N (x = 0,^ = 0) = 325 K within 
our saddle-point approximation. The inclusion of fluc- 
tuations away from the saddle-point can allow for the 
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more realistic value of J to be used. In the next section 
we discuss the effect of Sr doping and in the following 
we shall concentrate our studies on single crystals doped 
simultaneously with Sr and Zn. 



D. Effect of frustration - Sr doping 

In order to incorporate the effect of Sr doping in 
La2Cui_ z Zn z 04, we adopt the Shraiman and Siggia 
models In this model the holes introduced in the system 
(La2- 2; Sr 2 ;Cui_ 2 Zn z 04) via Sr doping are represented by 
an effective dipolar field that couples to the background 
magnetization current^ In Ref. [git was shown that this 
interaction leads to a reduction of the magnon gaps and 
spin stiffness, in agreement with the experiments^ As we 
shall now demonstrate, such reduction is not as strong if 
the compound is additionally doped with Zn, since in the 
presence of non-magnetic impurities the effective dipole- 
magnetization-current interaction should be multiplied 
by the bond dilution factor. This mechanism of reduction 
of frustration by nonmagnetic impurities has been con- 
sidered earlier by Korenblit et alA> within the isotropic 
0(3) NLSM. Here we discuss the role of anisotropies. 

In the Shraiman-Siggia model the interaction between 
magnons and the dipolar field in the absence of dilution 
can be written as 



Si, 



2A / dr 



d 2 rP, 



n x c^n, 



where 



•H.c, 



(13) 



(14) 



In the presence of dilution the above dipole-magnon 
interaction will also have to be modified following the 
procedure described in Appendix El Essentially, one can 
represent the semiclassical background spin distortion as 
a slowly varying SU(2) rotation of the Neel state R r . The 
hopping term for the doped holes involves the product 
ft-j- R,y, ^ and since these are rotations in neighbouring 
sites, i.e. R^R^", in the presence of dilution this product 
has to be replaced by 



R^R^ 



R,R, 



PiPj, 



(15) 



where i,j are nearest neighbor sites. Consequently, the 
coupling constant A between the dipolar field and the 
background magnetization current in Eq. Ill -it should be 
changed according to 



A -> KX, 



(16) 



because bl?>\) comes from 11150 . 12 Here K is the bond di- 
lution factor which in the homogeneous approximation is 

given by K =< p^pj >. 

A B n 



,'. x x % > 

' X <v 



clean system 



dirty system 



diffusive dipoles 



FIG. 3: Contribution to the magnon propagator from the 
dipolar fields. The figure on the left (A) contains the po- 
larization diagram of the doped holes. In (B) we introduce 
disorder through the scattering of the holes by impurities and 
this is represented effectively by the magnon-dipole bubble 
diagram (C) where we use the diffusive dipole propagator. 



P M is the dipolar field representing the spin current of the 
holes, \T/ is the spinor wave function of the doped holes 
with dispersion centered at (tt/2,±tt/2) and symmetry 
related points in the Brillouin zone, ~o are the three Pauli 
matrices, and /i is a lattice index. 



III. GENERAL MODEL 

The complete action describing the simultaneous effect 
of magnetic dilution by Zn and frustration by Sr reads 



J 



2gc(K/P a 



J At J A 2 r{{d T n) 2 + Z[c 2 (X7n) 2 +m 2 a n 2 a +m 2 c n 2 c ]}-2\K J At J d 2 rP, 



n x c^n + S d , 



(17) 



r 



where 



Sd 



-I At 



A 2 r P p G" 1 P M 



(18) 



describes the fluctuations of the dipolar field. Since the 
dipolar field is a pseudo-fermionic composite field, it is 
an operator that represents the spin current of the doped 
holes, the dipolar susceptibility (inverse dipole propaga- 
tor) can be obtained from a polarization diagram (see 
Fig.EJ. 



In what follows we adopt the same procedure used by 
Sachdev in Ref . [2a and we use the results from the Fermi 
liquid polarization diagram in two dimensions. However, 
while in Ref. [2a the ballistic limit (clean system) was used 
in the calculation of the polarization diagram, see left di- 
agram in Fig. El here we consider the diffusive limit (dirty 
system) because of the scattering of the doped holes by 
the Sr impurities, see middle diagram in Fig. 03 This as- 
sumption is consistent with recent results from Liischer 
et alm& who have shown that the dynamics of these dipo- 
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lar fields is highly diffusive. Thus, we shall write for the 
dipole propagator in the diffusive limit (see right diagram 
in Fig. EJ 



G D (q,uj n ) 



Dq 2 



Dq 2 + \uj n 



(19) 



where K d is proportional to the inverse static susceptibil- 
ity and D is the diffusion constant, which is assumed to 



be large. 



Because of the peculiar dipolar-magnon coupling from 
Eq. (I3J, the magnon propagator will be renormalized 
by the fluctuations of the dipolar field (see Fig. EJ) - The 
self-energy correction to the magnon propagator is 



E^(q,w n ) = K d (\ 



i^gcK 2 1 ^ f d 2 k 



(k + 2q) 2 



Dk 2 



> 2 w 2 + Zc 2 (k + q) 2 + Zm 2 Dk 2 + \u n \ 



(20) 



where e Q( 3 7 is the completely antisymmetric tensor. 

By summing up the one-loop corrections to the magnon Green's function, we obtain 

oo 

Gm = Gm J^(SmGm) j = [G A j — 1 - 

i=0 

Writing this expression explicitly, we get for the magnon propagator 

G%~\q, Lu n ) = (gcK/P^)- 1 ^ + Zc 2 q 2 + Zm 2 ) - S M (0) - 



q=0 



(21) 



(22) 



where we expanded the self-energy around zero momentum up to the second order term, in order to obtain a correction 
to the gaps and spin stiffness, and thus 



(gcK_ 



Kc ^9 2 S^(q,0) 
9 









q=0_ 


q 2 + 





Zm„ 



gcK 



£(o,o) 



2 dq^dqv 

From the above equation the expressions for the mass and the spin stiffness renormalizations are readily derived 

6^ 9 2 E^(q,0) 



M^(x,z) = Zm 2 a - ^-^(0,0) and p s (x,z)=Kp 8 



2 dq M dq„ 



q=0 



(23) 



(24) 



where c/g is the bare spin stiffness in a clean system. Explicitly, the mass renormalization due to Sr and Zn impurities 
reads 



where 



Ml{x,z) = Z [m 2 a - K d {gcXf^ZI{z)] , 



Dk 2 



2 w 2 + Zc 2 k 2 + Zm 2 Dk 2 + \w r , 



The spin stiffness, on the other hand, renormalizes according to the formula 

(k + 2q) 2 



(25) 



(26) 



(27) 



p s (x, z) = Kp s 



1 



2/3 



d 2 k / d 2 



(2tt) 2 I dq 2 u 2 + Zc 2 {k + q) 2 + Zm 2 



Dk 2 



q=0/ L>k 2 



The new momentum cut-off A for the dipoles, which renormalizes the spin stiffness, is set by Uf = \/ttx, because our 
theory should be valid at distances much larger than the average distance between the Sr impurities. 

In the zero temperature, T — ► 0, and highly diffusive, D — > oo, limits the expressions for M a and p s read 



M 2 a [x,z;T^0] = z\ml-K d {g C \) 2 e 2 Z^ 2 ^ 



Aire 4 



fc 2 c 2 \i/ 2 



(28) 
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and 



Ps[x, z;T -> 0] = ifp s < 1 - 



Kd ( 5 cA) 2 e2 /37 Z 1 /2 TO/3 



Aire 4 



1 *4c 

2 m% 



- 1 



3/2 



(29) 



IV. OBSERVABLES 



M(x, fixed z) ( ™ units ofj ) 



Eqs. l[28j) and (|29l are the most important results of 
this paper. As we can clearly see, although both M a 
and p s are reduced when only Sr (A ^ and Z — 1) or 
only Zn (A = and Z < 1) are doped independently into 
La2Cu04, when combined a nonmonotonic and reentrant 
behaviour can indeed occur due to the X 2 VZ coeficient 
in the self-energy correction (see Eqs. H28I1 and J2HJ)). 

In what follows we discuss the evolution with both Zn 
and Sr doping of several spectral, thermodynamic, and 
magnetic properties of the La2- a; Sr a ;Cui_ 2 Zn 2 04 system 
using the results obtained from the previous section. As 
we shall see, the nonmonotonic behaviour of different 
physical observables, such as the order parameter, the 
Neel temperature, and the weak-ferromagnetic moment, 
all follow from the nonmonotonicity induced by the com- 
petition between dilution and frustration in M a and p s . 



A. Dzyaloshinskii-Moriya or in-plane gap 

Juricic et al& have recently studied the evolution of 
the DM gap with Sr in I^-^Sr^CuC^ (see also Liischer 
et alm^Y. One of the most important findings reported 
in Ref. |6j was that the DM gap gives robustness to the 
canted Neel state and vanishes at a critical concentra- 
tion given by x = const. (Adm/ J) ~ 2%, where const. 
is O(l) (see also curves in Figs. Q]andQiJ|. The mecha- 
nism for the reduction of the DM gap was traced back to 
the self-energy corrections due to the dipolar frustration 
introduced by Sr doping, M a with Z = 1 (z = 0), and 
the theoretical curved was shown to agree quite well with 
the one-magnon Raman experiments of Gozar et alw& 

In fact, these experiments, performed at 10 K, show 
that at x = z = the Dzyaloshinskii-Moriya (or in- 
plane) gap is m a = 17.5 cm -1 (or 2.16 meV), whereas 
it is reduced by almost 30% at x = 1%, and it van- 
ishes at x = 2%. In Fig. we exhibit the Sr dependence 
of the DM gap for various fixed Zn concentrations (the 
DM gap is given in units of J). As it is evident from 
the plot, for small x and z the DM gap is always re- 
duced. At large x, however, a reentrant behaviour is 
observed already for very small Zn concentration. This 
effect is largest close to the Zn-free critical concentration 
x = 2%. In fact, while for x = 2% and z = the DM 
gap is zero^Afi, we find that for x = 2% and z = 15% the 
DM gap is A DM w 7.5 cm -1 , right above the lowest low- 
energy cutoff for one-magnon Raman spectroscopy (~ 5 
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FIG. 4: Dependence of the DM gap (Adm = M a , in units 
of J) on Sr doping x at various fixed Zn concentration: 
z = 0,0.05,0.1,0.15,0.2. Although Zn and Sr independently 
reduce M a , when combined an enhancement of M a does oc- 
cur. 
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FIG. 5: Dependence of the DM gap (Adm = M a , in units 
of J) on Zn doping z at various fixed Sr concentrations: x = 
0, 0.01, 0.0198, 0.0202, 0.03. The nonmonotonic and reentrant 
behaviours for M a is evident. 



cm -1 ). Thus, we predict that low-energy Raman exper- 
iments in La2- £! ;Sr a ;Cui_ z Zn z 04 samples with x w 2% 
and z — 15% should observe a clear signal of a one- 
magnon mode with energy close to 7.5 cm -1 in the Bi ff 
scattering geometry ; 11 ^ 7 such as the ones performed by 
Gozar et alm& in Zn-free samples. 
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Fig. shows the evolution of the DM gap as a func- 
tion of Zn doping for various fixed Sr concentrations. 
The nonmonotonic (for x = 0.0198) and reentrant (for 
x > 0.02) behaviours predicted by our theory for the 
DM gap are also evident from this plot. As expected, 
all curves collapse into a single curve in the highly Zn- 
diluted regime. For completeness we show in Fig. 0a 3D 
plot with the evolution of the DM gap with both Zn and 
Sr doping. 

As we have discussed in the previous section, the key 
mechanism for such nonmonotonic and reentrant be- 
haviour observed in the DM gap, through dilution by Zn, 
is the decrease of the effective dipolar-magnon coupling 
constant A — > KX in Eq. itTfi)) . which therefore reduces 
the self-energy corrections to the magnon gaps, see (Eq. 
El- 
Let us now list a number of predictions from our stud- 
ies, for a few selected doping concentrations of Sr and 
Zn, that can be verified experimentally. 

1. For La2Cuo.96Zno.04O4, that is x — and z — 0.04 
where \[Z ss 0.9, we estimate that the reduction of 
the DM gap will be of about 10% from the undoped 
x — z — value. 

2. For La1.99Sro.01Cuo.97Zno.03O4, that is x — 0.01 
and z = 0.03 where \[Z « 0.92, the reduction of 
the DM gap will be smaller than the naive 8%, due 
to Zn, over the 30% due to Sr. This is because 
of the dilution of frustration in the dipole-magnon 
coupling. 

3. For La1.9sSro.02Cuo.85Zno.15O4, that is x = 0.02 
and z = 0.15 where \/~Z ~ 0.6, we predic that the 
Dzyaloshinskii-Moriya gap is actually nonzero. In 
fact, due to the reduction of 40% in the dipole- 
magnon interaction, the DM gap is found to be 
7.5 cm -1 and thus large enough, so that it can be 
accessed with one-magnon Raman spectroscopy. 



M(X,0) (in units of J) 




0.005 



0.01 0.02 0.03 0.04 0.05 0.06 
Sr content (x) 

PIG. 7: Dependence of the XY gap (Axy = M c , in units of 
J) on Sr doping x at zero Zn concentration: 2 = 0. 



B. XY or out-of-plane gap 

The mechanism for the reduction of the XY gap is the 
same as the one discussed in the previous section, see Eq. 
II28II . For the XY gap the available data also comes from 
one-magnon Raman spectroscopy experiments, 1 ^ 1 For x = 
z = and at 10 K the XY gap is m c — 36 cm -1 (or 4.3 
meV), and is reduced by almost 15% at x = 1%. 

The dependence of the XY gap on Sr doping, according 
to our theory, is shown in Fig.0 We find that at x = 1% 
and z = the XY gap is reduced by almost 10%, which 
is not far from the experimentally measured value, and 
it vanishes at x = 0.04. When Zn is also included in 
the calculations we obtain a similar nonmonotonic and 
reentrant behaviour, as observed for the DM gap (see 
the 3D plot in Fig. EJ) . 

The vanishing of the XY gap at x = 4% for z = 
has another very interesting consequence. As it has been 
proposed by Juricic et al.fi for z = and x > 0.02, 



M^X,Z) < in units of J ) 



M(X.Z) < in units of J ) 
C 





Zn content (z) 




Zn content (z) 



Sr content (x) 



FIG. 6: 3D plot with the dependence of the DM gap (Adm 
M a , in units of J) on both Sr and Zn doping. 



FIG. 8: 3D plot with the dependence of the XY gap (Axy 
M e , in units of J) on both Sr and Zn doping. 
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after the DM gap has disappeared, the magnetism be- 
comes incommensurate (the staggered moment becomes 
helicoidal along the 6-axis) with an incommensurate wave 
vector Q || 6 and with magnitude 

Q = yJx*-M c (x,0y, (30) 

where M c = M c /p s is dimensionless. As a consequence 
of the nonzero character of M c , the incommensurability, 
for 0.02 < x < 0.04, deviates from the linear behaviour 
Q — x, as indeed observed experimentally^ For higher 
Sr doping, however, when M c {x > 0.04) = 0, the linear- 
relation Q = x becomes exacts 

C. Order parameter 

Because of the DM and XY anisotropies, at zero ap- 
plied magnetic field the staggered order parameter is 
oriented along the orthorhombic b axis with magnitude 
(Xq(x, z). In the presence of a longitudinal magnetic field, 
B || b, the spins start to rotate on the be plane as de- 
scribed theoretically by Silva Neto and Benfatto in Ref. 
Ull (see also Benfatto et al. in Ref. H% and experimen- 
tally measured with neutron diffraction by Reehuis et al. 
in Ref. At a critical field determined by the DM gap 

H c (x,z)=M a (x,z), (31) 

a spin-flop transition occurs, where the in-plane compo- 
nent of the order parameter becomes oriented along the 
orthorhombic a axis. Thus, any nonmonotonic and reen- 
trant behaviour of the DM gap will result on a similar 
reentrant behaviour for the critical field. This is a pre- 
diction that could be detected by neutron scattering ex- 
periments. 

D. Neel temperature 

In the case of single crystals doped with Sr and Zn, the 
available data from Ref. were obtained with samples 
which contained less than 30% of Zn. Therefore, in this 
doping range the dependence of Tn(z) is linear and is well 
described with the bond percolation factor K{z) — 1— 3z. 

The Neel temperature, Tff{x,z), where the order pa- 
rameter vanishes, is given by Eq. illlfl . with the stiff- 
ness and magnon gaps renormalized by the impurities, 
see Eqs. J2EJ and ffity. As it happens with the DM gap, 
we find that T/v exhibits a nonmonotonic and reentrant 
behaviour as the Zn concentration is increased, see Fig. El 
In particular, we find a monotonic decrease in the slope of 
the curves Tn(x,z = 0) and Tn{x,z = 0.15), in qualita- 
tive agreement with the experiments of Hiicker et alA, see 
Fig. Since we neglect the self-consistent renormaliza- 
tion of the magnon gaps with temperature, the agreement 
of our theoretical curves for T N with experiments is only 
qualitative, in contrast to the low-temperature results for 
the DM gap discussed above, which are quantitative. 
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FIG. 9: Dependence of Tn (in K) on Zn doping z at various 
fixed Sr concentrations x = 0,0.01,0.0198,0.0202,0.03. The 
nonmonotonic and reentrant behaviours are evident. 
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PIG. 10: Dependence of Tn (in K) on Sr doping x at two 
different Zn concentrations: z = and z = 0.15. The slope 
of the curves decreases monotonically with increasing Zn con- 
tent. 



It is worth emphasizing that the nonmonotonic be- 
haviour exhibited by our theoretical curve x = 0.0198 in 
Fig. is experimentally observed already at x = 0.017. 
Moreover, the 2D-Ising-like behaviour of the curves in 
Fig. E| are actually an artifact of the dimensionality (we 
are considering an easy-axis 2D NLSM) and of the ap- 
proximation (we are neglecting the thermal renormaliza- 
tions of the gaps) . We expect that by including the self- 
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consistent thermal renormalizations and interlayer cou- 
pling, the agreement between theory and experiment in 
Fig- El will be satisfactory. Nevertheless, the qualitative 
agreement with experiments gives strong support to the 
model and mechanisms considered here. 
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E. The weak-ferromagnetic moment 



APPENDIX A: DERIVATION OF SITE AND 
BOND PERCOLATION FACTORS 



The staggered pattern of tilted octahedra in 
La2Cu04 is known to be responsible for a weak- 
ferromagnetism, signaled by a cusp in the low-field mag- 
netic susceptibility^ Within the language of the NLSM 
(see Ref.|§), such weak-ferromagnetic moment is propor- 
tional to the Neel order parameter via the equation 



1 

2J 



(n) x D+] 



(32) 



Since within the Neel phase the average value of the stag- 
gered magnetization is (n) = (0, <tq,0) (in the abc or- 
thorhombic coordinate system), any nonmonotonic and 
reentrant behaviour observed in ao will cause also a sim- 
ilar effect in the weak-ferromagnetic moment, and this 
can be accessed in magnetic susceptibility experiments. 



V. CONCLUSIONS AND OUTLOOK 

In this paper we revisited the problem of the dilution of 
frustration in La2- a: Sr a ;Cui_ z Zn 2: 04, within the frame- 
work of a generalized NLSM that includes DM and XY 
anisotropics. We showed that dilution by Zn weakens 
the frustration by Sr through the reduction of the dipole- 
magnon coupling constant, see Eq. iflfi)) . This leads to a 
nonmonotonic and reentrant behavior not only for T/v 
but also for other observables like the order parame- 
ter, the weak-ferromagnetic moment, and the anisotropy 
gaps. 

Most remarkably, we predict that for x ~ 2% and 
z = 15% the DM gap is approximately 7.5 cm -1 , that 
is, larger than the lowest low-frequency cutoff for Raman 
spectroscopy (~ 5 cm -1 ) and thus likely to be observed in 
one-magnon Raman scattering. Furthermore, when the 
WL expression for the bond percolation factor is incorpo- 
rated into our NLSM description, not only as a reduction 
factor for the spin-stiffness but also and most importantly 
for the reduction of the anisotropy gaps, we find that 
our NLSM with dilution describes correctly the data for 
T/v (x = 0, z), also in the highly Zn-diluted regime. Fi- 
nally, we have also found that the XY gap vanishes, in 
the absence of dilution, for x = 0.04 and this is consistent 
with the deviation from linearity, for 0.02 < x < 0.04, 
of the incommensurate peaks seen in neutron scattering 
within the spin-glass phase of La2- a: Sr a: Cu04. 



This derivation can be found also in Ref. |23, but for 
the sake of completeness we will here go through the basic 
steps of the derivation. 

Averaging the site percolation factor we obtain 



Poo Poo ( % ) — 

= 1 



1 



Na 2 
Nzn/N =l-z, 



d 2 rp(v) = 



a 2 {N-N Zn ) 



Na 2 



(Al) 



where averaging p(r) over a 2D volume implies that we 
have to normalize it with Na 2 , N being the number of 
sites in the Cu-lattice. 

We consider pj to be smooth, and thus we can expand 
it in the neighborhood of the z-th site pj — pi + aVp./. 
Hence, the bond percolation factor Kij = PiPj in the 
continuum limit becomes 

K(r) = p(r)[p(r) + aV>(r)] = p{r) + -aV>(r), 



leading to 



K 



K{z) 

1-2- 

1- z- 



1 



Na 2 
a 

2Na 2 
a 

2Na 2 



p(r) + -aVp(r) 



d 2 r 
d 2 rVp(r) 



(— 1)V^ / contours 

Zn 



a7V z „4a 

= 1 — z — — „ = 1 — iz 

2Na 2 



(A2) 



Using Stoke's theorem in the above formula we obtained 
the integration over a contour of the 2D volume, which 
splits into the sum of contours over Zn impurities. The 
minus sign accounts for the opposite orientation of the 
small contours with respect to the larger one. 



APPENDIX B: NLSM WITH DILUTION 

In La2CuC"4 the DM vectors are in good approxima- 
tion perpendicular to the Cu-Cu bonds and change sign 
from one bond to another, while the XY matrices pro- 
vide an easy-plane anisotropy. It is worth noting that 
the pure 2D system defined by the action (QJ does not 
display a rotational symmetry, so it can have order at fi- 
nite temperature without violating the Mermin- Wagner 
theorem. 
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Before we proceed with the derivation of the NLSM, 
let us recall some technical details: d + is a vector in the a 
orthorhombic direction; and i is an index on a 2D lattice, 
so i should be understood as (p, q) and (— l) 1 as (— l) p+q . 
Furthermore, we shall use that 



E D y = E D ^+ D ^ = 2 E d 



+,») 



(Bl) 



<hs> 



a = (-l) 4 n 4 + ah - -(-ly^ndl (B2) 



with nj • Ij = 0, 



n 3 = n, - reditu + -r^d^nu (B3) 



and finally 



^-1 



d 2 r. 



(B4) 



Let us first transform the Heisenberg term of the 
Hamiltonian. Using Eqs. <IB2f) and <IB3f) and neglecting 
0(a 3 ) terms we get 



H H = JS 2 E a i n 

<i,j> 

= ^ 2 E 

<i,j> 

= JS 2 a 2 E 



(Vn.^ + 41| 



(B5) 



In the continuum limit the Heisenberg Hamiltonian reads 



JS 2 



d 2 r [(Vn) 2 + 81 2 ] 



(B6) 



In the diluted case, we should multiply fij with pi. 
Hence, Eq. IB 51 would be modified as follows 

H d H = JS 2 E 

<i,3> 

where the index "d" stands for "diluted". We want to 
treat the static impurities as an average effect, thus we 
substitute piPj in the above equation with < pipj >= 
K(r), the bond percolation factor (for more details, see 
chapter 3 or Ref . 23) . We simplify the problem even more 
by considering the averaged one, taking < K{v) >= K(z) 
(for the shortening of the notations in what follows we 
will use K instead of K(z)), where z is the Zn concentra- 
tion. Thus, in the diluted case Eq. IIBfill reads 



Following a similar procedure, we may transform the DM 
and XY terms of the Hamiltonian. For the DM Hamil- 
tonian we get 



TTCl 

11 DM — 



45 ^|d 2 r[d + .(nxl)]. 



(B8) 



For the XY Hamiltonian in the continuum limit we find 



H XY — 



2S 2 K 



J d 2 r [(Ti - T 3 )nl] , (B9) 



where we neglected the small terms like r(Vn) 2 and Tl 2 . 

Now, we will discuss in detail a Wess-Zumino term 
in 2D for the diluted case, since it didn't appear in the 
literature. For the ID case we refer the reader to Ref. 
(clean system) and Ref. I23I (diluted system). 

Following Fradkin (see Appendix A or Ref. we 
write the Wess-Zumino action on a lattice (notice that 
(p, q) are the indices along x and y directions respec- 
tively, in contrast with i and j, which take values on a 
2D lattice) 

sswz = 5n-nxd Q n (bio) 

S / dx ^2Swz[fl(p, q)} 

■'" P,9 

V*/2 N y 

E E {W2p, <?)] + s[n(2 P - 1, q )]} 



Swz 



S 



p=l g=l 



N x N v /2 

+ E E ^ n (p> 2 ^ + 5 ["(p> 2 « - !)]> 

p=l q=l 



(BH) 



In the second and third lines of Eq. ijBllll we recognize 
SSwz, which can be expressed through CI. Using the 
spin decomposition (IB2I) . we get 

6 x Sl(p,q) = n(2p,g) + n(2p-l,g) 

= (-l)*»[n(2 Pi q)-n(2p-l,q)] 

+ 2al[2p, q) + 0(a 2 ) 

= ad x n + 2a\. (B12) 

Analogously, we transform the third line of Eq. ijBllll 

S y n(p, q) = ad y n + 2al. (B13) 

Let us note that 

Nx/2 Ny 

I . — * I . — * I / 

(B14) 



— 9 — 9 



p=l q=l 



all sites 



1 

2^ 



d 2 r. 



Plugging JbToIi , CHS} and l|BT3)l into ifBTTT and hav- 
ing in mind l|B14l) . we obtain for the Wess-Zumino action 
(decomposing ijBlflll into staggered and uniform compo- 
nents, we keep only first order terms in a) 



K = 



JS 2 K 



J d 2 r [(Vn) 2 + 81 2 ] . (B7) Swz = T J dX ° E[( Vn < ; + 41 ^ ' ( n * x d o n i)l ( B15 ) 
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The first term in Eq. l|B15ll is a topological term. It was 
demonstrated by HaldaneSi that in D>1 this term sums 
to zero in the AF background. 
The total Euclidean action reads 

Se = —iSwz + Sh + Sdm + Sxy- (B16) 
Now we are ready to write a Wess-Zumino Lagrangian 

Ly/z 



J d 2 r[l-nx d T n}. (B17) 



In the presence of dilution the action l|B17J) will be mul- 
tiplied with Poo. The explanation is the following: in 
the expression for SSwz (IB 10p we will have a factor p 3 , 
which by definition is equal to p (and we further simplify 
the problem by taking < p >= Poo). We can carefully do 
the procedure in Eq. IIBllj) in the diluted case, but since 
we neglect 0{a 2 ) terms, the answer will be the same. 

Let us now summarize the results obtained so far. The 
total Euclidean action in the diluted system reads 



S = drL t 



(B18) 



with L to t = L W z + L H + L DM + L X y and 



Lwz — 
Lh = 
Ldm 



■ J d 2 r [1 • (n x 9 T n)] , (B19) 



JS 2 K 



J d 2 r [(Vn) 2 + 81 2 ] , 
AS 2 K f , 2 

— — / d 2 r [d+ ■ (n x 1)] , 



L 



2S 2 K 



XY 



|d 2 r[( ri -r 3 )n 2 ]. 



Now we integrate out 1 in a sense of a saddle point solu- 
tion. We have to find a solution of an equation 



51 



and plug it into Eq. (IB18J1 . Doing this, we get 



BJSaK 



(n x d T n) 



2Ja 



(n x d 



+), 



(B20) 



and 



S 



2gcJ J [K 
KB 2 + n 2 a +KT c n 2 c ], 



(d T n) 2 +Kc 2 (Vnf 



(B21) 



where we defined 



gc = 8Ja 2 - (bare) inverse transverse susceptibility, 
c = 2V2JSa - (bare) spin-wave velocity, 
D+ = ^/2gcS 2 /Ja 2 d+ = 2V2Sd~e a - DM vector, 
T c = {AgcS 2 /a 2 ){T l - T 3 ) - XY anisotropy. 

It is convenient to introduce the notations: D + = m a and 
r c = m 2 . Plugging this into Eq. JB211 and rewriting it 
in a conventional way we get a final expression for the 
total action in the presence of dilution 



S = 



2gcK/P a 



■ J dr J d 2 T [(d T nf + Z{c 2 (Vn) 2 + m 2 a n 2 a + m 2 n 2 }] , 



(B22) 



where Z = K 2 /P„ 
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